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Abstract. We investigate the bright cut-off of the [O m] /15007 planetary nebula luminosity function (PNLF), that has been sug- 
gested as a powerful extragalactic distance indicator in virtue of its observed invariance against populations effects. Theoretical 
PNLFs are constructed via Monte-Carlo simulations of populations of PNe, whose individual properties are described with the 
aid of recent PN synthetic models (Marigo et al. 2001), coupled to a detailed photoionisation code (CLOUDY). The basic depen- 
dences of the cut-off magnitude M* are then discussed. We find that : (i) In galaxies with recent or ongoing star formation, the 
£ — | modelled PNLF present M* values between -4 and -5, depending on model details. These are very close to the observationally- 

(f) , calibrated value for the LMC. (ii) In these galaxies, the PNLF cut-off is produced by PNe with progenitor masses of about 2.5 

M e , while less massive stars give origin to fainter PNe. As a consequence M* is expected to depend strongly on the age of the 
last burst of star formation, dimming by as much as 5 magnitudes as we go from young to 10-Gyr old populations, (iii) Rather 
than on the initial metallicity of a stellar population, M* depends on the actual [O/H] of the observed PNe, a quantity that may 
differ significantly from the initial value (due to dredge-up episodes), especially in young and intermediate-age PN populations, 
(iv) Also the transition time from the end of AGB to the PN phase, and the nuclear-burning properties (i.e. H- or He-burning) 
of the central stars introduce non-negligible effects on M* . The strongest indication derived from the present calculations is a 
serious difficulty to explain the age-invariance of the cut-off brightness over an extended interval, say from 1 to 13 Gyr, that 
observations of PNLFs in galaxies of late-to-early type seem to suggest. We discuss the implications of our findings, also in 
relation to other interpretative pictures proposed in the past literature. 
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?H 1 1. Introduction In practice, the PNLF as a candle works in the following 

way: Once the flux in the [O in] /15007 emission line, F5007 
The Planetary Nebula luminosity function (PNLF) constitutes (i n cgs un i ts ), is measured for a sample of PNe in an external 
one of the most attractive indicators to set out the extragalactic galaxy, it is converted to m 50 07 magnitude, 
distance scale. From the observational side, PNe can be identi- 
fied in galaxies from the LMC to as far as the Virgo cluster, OT 5007 = -2.5 logF 500 7 - 13.74 (1) 
and are numerous in some of the key galaxies that link the (Jacoby 1989) The number of objects at any given ^ bin 
"local" distance scale given by Cepheids, with the more dis- constitutes the PNL F 5 jV( m5007 ), which can be fitted over the 
tant scale given by the surface brightness fluctuations (SBF) brightest ma g nitu des by the function 
method, the Tully-Fisher relation, and SNe la. The PNLF is 

suitable for determining distances as far as 20-25 Mpc, beyond N(m 5m ) = e 03Q7m5001 [l - e 3(m * "'" 5007) ] (2) 



which exposure times for reasonable S/N-ratios become pro- 
hibitively long. The importance of the PNLF for the cosmolog- 
ical distance ladder is fully described by e.g. Ciardullo (2003), 
Ciardullo & Jacoby (1999) and Jacoby, Ciardullo & Feldmeier 
(1999), Jacoby (1997), to whom we refer for all details. 



(see Ciardullo et al. 1989). The key quantity is the cut-off mag- 
nitude m*. According to several arguments (see below), its cor- 
responding absolute value, M*, appears to be constant in most 
galaxies, hence providing an excellent (secondary) standard 
candle. Therefore, a measure of m* gives a distance estimate 

in a single step, provided that interstellar reddening is suitably 

Send offprint requests to: Paola Marigo taken into account. The present-day calibration of AT, based 

e-mail: marigo@pd.astro.it on a sample of galaxies with SBF and Cepheid distances, and 
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correcting for a small metallicity dependence, is M* — -4.47 
(Ciardullo et al. 2002). 

The possibility of systematic variations in the value of M* 
among galaxies has been investigated in several different ways. 
From the observational side, the evidences for a nearly con- 
stant maximum luminosity for PNe are numerous, as reviewed 
by e.g. Ciardullo (2003) and Jacoby (1997). The main points 
are: (1) m* does not significantly change with galactocentric 
radius in the galaxies M31, M81, NGC4494 and NGC5128, 
where significant population gradients should be present. (2) 
There is no hint of large systematic variations in m* among 
galaxies of different types in the galaxy clusters/groups of M8 1 , 
NGC 1023, NGC 5128, Fornax, Leo I, and Virgo. Overall, this 
is indicating that age and metallicity variations have little effect 
on M*. On the other hand, comparison with Cepheid distances 
reveals a small non-monotonic dependence on metallicity, ex- 
pressed by the system's oxygen abundance, which can be com- 
pensated for by a relation quadratic in [O/H] (Ciardullo et al. 
2002). 

Although these observational arguments seem compelling, 
they are not accompanied by a solid, unique theoretical frame- 
work. In other words, different authors do not agree on the 
reasons why M* should be constant, reaching seemingly con- 
tradictory conclusions over the years. First, early theoretical 
arguments by Jacoby (1989) and PN modelling by Dopita et 
al. (1992) explain why the [Om] /15007 flux depends little on 
the PN metallicity. At the same time, the PNLF simulations 
by Mendez et al. (1993) indicate that a substantial dependence 
on population age might be present, as indicated by the dif- 
ferences in the best fits of the LMC and Milky Way PNLFs. 
Dopita et al. (1992) also suggest that a dependence on age 
should be present. While Mendez et al. (1993) claim that the 
bright end of the PNLF is populated by optically-thin objects, 
Dopita et al. (1992) model optically-thick objects only and, af- 
ter metallicity and population corrections, succeed in reproduc- 
ing the PNLF and the constant cutoff brightness. Similar open 
questions concern the importance of hydrogen- (Mendez et al. 
1993) vs. helium-burning (Ciardullo et al. 1989; Jacoby 1989) 
central stars. 

More recently, the apparent lack of predicted population 
age effects is explained by Ciardullo & Jacoby (1999) by com- 
pensating circumstellar extinction. Young and massive central 
stars should be very luminous and result in over-luminous PNe 
(i.e. beyond the standard M*), which is not observed. The pro- 
posed solution by Ciardullo & Jacoby (1999) is that, since 
the most massive asymptotic giant branch (AGB) stars likely 
become dust-enshrouded because of heavy mass loss, the re- 
lated circumstellar extinction would reduce the [O m] /15007 
emission from these objects, always to at least M* or below. 
Provided that PNe with lower-mass progenitors are present, 
then the observed cut-off magnitude would be M*. If this is the 
case, one has yet to understand which is the typical age of the 
PNe that contribute to the cut-off, and explore the metallicity 
effects in the relevant age intervals. 

The theoretical investigations of the PNLF, quoted above, 
have different levels of sophistication in the models. Because 
of the complexity of the problem, simplifications and severe 
assumptions about the properties of PNe and their progeni- 



tor stars have always been introduced, like for instance: an 
initial-final mass relation (hereinafter also IFMR) independent 
of metallicity (e.g. Mendez et al. 1993), either optically-thin or 
thick PNe (Dopita et al. 1992; Mendez et al. 1993; Stanghellini 
1995; Mendez & Soffner 1997), a narrow mass range being re- 
sponsible for the bright end of the PNLF (Jacoby 1989; Dopita 
et al. 1992), and the restriction to either H- or He-burner central 
stars only. All these assumptions can be improved in some way, 
by using the full results of stellar evolution and the population 
synthesis theory. 

In the present paper, we elaborate on the theoretical frame- 
work for the PNLF predictions, based on our own simulations. 
The starting point for such a project has been the construc- 
tion of an updated synthetic code, that follows the PN evo- 
lution as a function of basic stellar parameters, such as pro- 
genitor mass, metallicity, type of post- AGB evolution, etc. The 
complete code is described in Marigo et al. (2001; hereafter 
Paper I), together with basic comparisons with the available 
data for Galactic PNe. Note that this code is much more com- 
plex than those previously used to investigate the PNLF, as it 
takes into account the chemical composition and wind velocity 
of the ejected material as predicted by AGB models, the inter- 
action of the multiple winds as the PN expands, the effect of 
shell thickening due to ionisation, etc. Moreover, there is no ar- 
tificial assumption about the conditions of optical thickness of 
the nebulae. The underlying model used to produce theoretical 
PNLFs is therefore more realistic and much less dependent on 
explicit or hidden assumptions than previous approaches. We 
will recapitulate the key features of our theoretical model in 
SectionlTTI 

Despite all improvements, the Marigo et al. (2001) mod- 
els still treated some characteristics of the ionized nebulae in 
a simplified way. For instance, the electronic temperature T e 
was introduced as a parameter and kept constant during the 
complete PN evolution. The emission line fluxes were then de- 
rived by simple formulas based on this temperature, on ele- 
ment abundances, and on the estimated sizes of the Stromgren 
spheres for different ions. Since the flux of the optical for- 
bidden lines - and in particular the [O ra] /15007 one - de- 
pends strongly on T e and the actual amount of ionized mate- 
rial, this simplification is clearly inadequate for an investiga- 
tion of the PNLF. Hence, in this paper (Section [2. 21 we relax 
the above-mentioned approximations by incorporating the pho- 
toionisation code CLOUDY (Ferland 2001) into our PN mod- 
els, thereby improving them further, and particularly towards 
the direction of accurate predictions of spectral features. 

Then, in Sect. [3] we describe the generation of simulated 
PN samples in galaxies, firstly performing some basic compar- 
isons with the empirical data of Galactic PNe. We show how 
present models can reproduce, to a large extent, some well- 
known properties of observed PNe, such as the correlations be- 
tween the ionized mass and nebular density with the radius, and 
those between particular emission line ratios. 

In Sect.@]we move to analyse the predicted synthetic PNLF, 
with the aim to discuss the expected behaviour of the bright cut- 
off, M*. We will in particular focus on the nature of the objects 
responsible for the bright cut-off, analysing the dependence of 
M* on physical and technical factors. 
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Finally, Sects. [5] and [6] summarise the results of present 
models in terms of basic dependences of the cut-off, and draw 
a few conclusions about the PNLF method, at the same time 
highlighting the main points of discrepancy between our find- 
ings and the picture used so far to interpret the PNLF results. 

2. Synthetic model for PN populations 

2.1. Overview of the PN model 

We follow the evolution of a single PN by using the synthetic 
code from Paper I. The main aspects of our PN models are: 

- The stellar evolution prior to the PN-phase is taken for a 
large grid of initial masses M; and metallicities Z from 
Girardi et al. (2000) for phases up to the first thermal pulse, 
and from Marigo et al. (1999) and Marigo (2001) during 
the complete TP-AGB phase (the AGB-phase with ther- 
mal pulses occuring). The TP-AGB evolution includes the 
main phase of mass loss, which later gives origin to the PN 
shell. We use the model predictions to follow the evolution 
in time and radius of the wind ejecta and their chemical 
composition. 

- We use post-AGB evolutionary tracks from Vassiliadis & 
Wood (1994) to model the evolution of the PN central 
star (CSPN). Such tracks are available for both H- or He- 
burning regimes. A transition time, t tr , is assumed between 
the end of the TP-AGB and the beginning of the post-AGB 
track at T eB = 10000 K. The CSPN tracks allow the com- 
putation of quantities such as the flux of UV ionising pho- 
tons and the velocity of the fast wind. 

- The PN shell is ionized by the UV photons and expands 
due to its initial ejection velocity. As a major improvement 
of synthetic models we include the following dynamical ef- 
fects: the shell shaping by the inner fast wind emitted from 
the hot CSPN (following Volk & Kwok 1985), the succes- 
sive engulfment of external layers expelled at early times 
with lower expansion velocities, and the interaction of mul- 
tiple winds, including the effect of shell-thickening due to 
ionisation (based on Kahn 1983). 

In the original model (Paper I), the description of ionisa- 
tion is simplified. The electron temperature T e in the nebula is 
assumed to be constant and close to 10 000 K for the entire H + 
region. Furthermore, the sizes of the H + , He + , He ++ regions 
are derived from simple considerations about the sizes of the 
different Stromgren radii, i.e. the radii below which the corre- 
sponding ionising photons would be responsible for complete 
ionisation. Moreover, the size of the ++ region, whence the 
[O m] A5QQ1 line is emitted, is simply assumed to be equal to 
the He + region; this latter assumption is based on the similar- 
ity between the He + and ++ ionisation potentials (54.4 and 
54.9 eV, respectively). Of course, these assumptions directly 
affect the predictions for emission lines such as He n /14686 and 
[Om] /15007. Despite these simplifications, it is demonstrated 
in Paper I that our PN models can reproduce quite well the ob- 
served correlation between the line ratios He n /14686/HyS and 
[Om]^5007/HeiU4686. 
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Fig. 1. The ionisation fractions of H, He, and O ions, as a func- 
tion of radius inside the ionised PN shell. The ionisation struc- 
ture has been computed with the aid of CLOUDY, as explained 
in Sect. 12.21 This particular PN model corresponds to the case 
M = 2.0 M@, Z = 0.008, t b = 1 000 yr, H-burning track, at an 
age of t — 3 000 yr, which is close to its maximum [O ra] A50QTI 
luminosity. Other relevant quantities are Mcspn = 0.685 M Q , 
M ion = 0.091 M , Rm = 0.059 pc, logr eff = 5.218, and 
log L/L e = 3.867. 



Nonetheless, these assumptions clearly are not adequate for 
an accurate modelling of the absolute [O in] /15007 flux from 
PNe, which is our main target in the present paper because: (1) 
the flux of any collisionally-excited forbidden line is very sen- 
sitive to the electron temperature T e , that therefore should not 
be assumed constant, but rather it needs to be consistently com- 
puted from the equations of energy balance inside the nebular 
material; and (2) according to detailed photoionisation models, 
the ++ ionised region is neither expected to have sharp edges, 
nor is the ++ fraction expected to follow so closely the He + 
one. These points are illustrated in Fig.^ 

To overcome the limitations of Paper I in this respect, the 
proper choice is to apply a complete photoionisation code to 
our evolutionary models. This step is described in the following 
section. 

2.2. Incorporation of CLOUDY 

We opt for the photoionisation code CLOUDY (Ferland 2000, 
2001), version 94, which is a publicly-available, widely-used 
and well-tested tool for describing the structure of ionised re- 
gions and their emission spectra. CLOUDY and other similar 
photoionisation codes have already been applied in the evolu- 
tionary modelling of PNe by a number of authors, among them 
Volk (1992), Stasiriska (1989), Stasinska & Szczerba (2001) 
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Dopita & Meatheringham (1990), Dopita et al. (1992), Dopita 
& Sutherland (2000), Phillips (2000). However, in all those 
cases the treatment of the underlying stellar population and of 
the nebula expansion was very crude, compared to our method 
of Paper I. 

In order to couple CLOUDY to our evolutionary code, we 
proceed in the following way: 

1 . The former synthetic code (Paper I) is run for a guess elec- 
tron temperature of T° = 10 000 K. 

2. At each time step, the radial structure (density, velocity, 
chemical composition, etc.) of the ionised shell is taken 
from the synthetic code and provided as input to CLOUDY. 

3. CLOUDY is then used to solve for the energy balance, ion- 
isation structure, and emissivities as a function of radius. 
From the CLOUDY output, we extract the most relevant 
information such as the volume-averaged electron temper- 
ature T e and the outcoming fluxes in some of the main emis- 
sion lines. 

To simulate the nebular structure given to CLOUDY a num- 
ber of parameters must be specified. For the present work, 
we adopt an expanding-sphere geometry with covering fac- 
tor equal to 1 (see Ferland 2001 for details), inner and outer 
radii R$ and 7?; on , and uniform hydrogen density N(H) (refer to 
Paper I for the precise definitions of the quantities). All the 30 
chemical elements allowed by CLOUDY are considered. For 
most elements, a scaled-solar abundance, corresponding to the 
adopted Z value, is assumed. The adopted fractions of H, He, 
C, N, O, Ne and Mg, on the contrary, are the ones predicted 
from the evolution of the PN progenitor, with all dredge-up 
events as well as hot-bottom burning being followed in de- 
tail (see Marigo 2001). It should be remarked that the abun- 
dances of these elements in the main expanding PN shell may 
change with time, due to the accretion of successive outer lay- 
ers ejected from the progenitor star at previous stages of its 
nucleosynthetic evolution. Note that such a variable PN com- 
position is being considered for the first time in a theoretical 
model for the PNLF. 

As for the central star, L and T e g values are converted into 
photon fluxes by using the following libraries of intrinsic stellar 
spectra, already available in CLOUDY: ATLAS9 from Kurucz 
(1993) for T eB < 50 000 K, and the PN nuclei non-LTE models 
from Rauch (1997) for T eB > 50 000 K. 

2.3. Consistency and accuracy checks 

Initially, our intention was to iterate over the sequence of oper- 
ationsnto|3]outlined at the beginning of Sect. 12.21 i.e. provide 
the output r e from CLOUDY as a new T® guess to the synthetic 
PN model, and iterate steps ^to|3]until convergence to a single 
value of T e had been reached. In practice, it has become clear 
that this iteration is not necessary, for the following reasons: 

First, the output T e values from CLOUDY turn out to be, 
in general, of the same order as the input value of 10000 
K, especially during the phase in which the highest HB and 
[O m] /15007 luminosities are attained. Typical values are com- 
prised between 8 000 and a maximum value of 15 000 K. Lower 
T e values are present only in the very initial phases of nebular 
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Fig. 2. Top panel: The volume-averaged electron temperature, 
T e , resulting from CLOUDY, as a function of time, for a se- 
quence of PN models computed for various values of the input 
guess T°, namely 8 000, 10000, and 12000 K (dotted, dashed, 
and continuous lines, respectively). The models correspond to 
the case M ; = 2.0 M @ , Z = 0.008, f tr = 1 000 yr, H-burning 
tracks. Bottom panel: The resulting absolute magnitude in the 
[O in] /15007 line, M5007, as a function of time, for the same PN 
models. M5007 is only little affected by the value of the input 
guess T%. 

ionisation (as the central star is heating up), and in very lat- 
est phases of recombination (as the central star fades on the 
white dwarf track). These stages of low T e present either very 
low nebular emissivities, or very low surface brightnesses, and 
likely do not correspond to observable PNe. 

Second, the input T l J provided to our synthetic code has a 
very small impact on the structure and evolution of the neb- 
ula, and consequently hardly affects the output T e values from 
CLOUDY. This is illustrated in Fig.|2]which presents the output 
T e from CLOUDY for a few models, corresponding to the same 
PN progenitor star but adopting three different values the input 
temperature T°, i.e. 8 000, 10 000, and and 12 000 K. Evidently, 
the output T e values differ just little between the various cases. 
For the main phase of PN evolution, in the age interval between 
2 000 and 4000 yr (where the [Om] /15007 luminosity attains 
the highest values), the T e values differ by less than 4%. The 
differences are caused by the modest changes in the density - 
and hence radial thickness - of the nebular shell, for which the 
input T® plays a minor role. 

For this same set of models, Fig.[2]also presents the result- 
ing absolute magnitude in the [Om] /15007 1 , M5007, as a func- 
tion of time. Again, the trial input 7* has just a modest impact 
on the line flux calculated by CLOUDY. The largest differences 

1 M 5007 is computed from Eq. Q, assuming a PN distance of 10 pc. 
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are of about 0.3 mag for the models illustrated in Fig.|5J we find 
that at larger stellar masses, that produce the brightest PNe, the 
differences are even lower. From tests we have performed, we 
expect that on average, by using a constant T® of 10 000 K, the 
typical PN models will have M5007 magnitudes as accurate as 
to about 0.15 mag. 

Such an accuracy may not seem good enough for the study 
of a potential standard candle. However, we should recall that 
this theoretical indetermination of M* is comparable to the 
global uncertainty affecting the observational derivation of m* 
(see Jacoby et al. 1990). Moreover, we expect that the errors 
in our predicted M5007 will be mostly systematic, whereas our 
study intends to explore mainly the differential effects - i.e. the 
changes in M* between different galaxies. We also remark that, 
although the absolute values of M* is of high interest, an accu- 
rate theoretical calibration of its value is still out of reach, as 
we will demonstrate in this paper. 

2.4. Summary of parameters 

The low sensitivity of our results to the input guess T®, effec- 
tively cancels this quantity from our list of input parameters. 
The adjustable parameters that remain in our models are: 

- the initial mass M\ and metallicity Z; 

- the transition time f tr ; 

- the choice of either H- or He-burning post-AGB track. 

Other important quantities like the initial He content, and the 
CSPN mass, are defined by the relationships already incorpo- 
rated in the synthetic evolution models, such as the helium-to- 
metal enrichment ratio (see Girardi et al. 2000), and the initial- 
final mass relation (see Marigo 2001), which has a dependence 
on metallicity. 

The typical values for t u , and the type of post-AGB evo- 
lution, have a great impact on the evolution of single PN. For 
a given mass and metallicity, they completely settle the evolu- 
tionary timescale of the central star, and hence the duration of 
the PN phase. It is worth recalling the present uncertainties in 
these key ingredients of our models: 

The transition time f tl is very badly predicted by theory. Its 
value depends on several details of the final phase of mass loss 
after the termination of the AGB-phase, when the star is get- 
ting hotter. Main factors driving f tr are the velocity by which 
the residual stellar envelope is peeled off from the star, and the 
velocity by which the same envelope is burnt from below by nu- 
clear reactions. Educated guesses about these processes can be 
made (see e.g. the discussion in Stanghellini & Renzini 2001), 
but lead to a general picture of great uncertainty, in which the 
only firm conclusions are that f tr should generally decrease with 
stellar mass, and have values of the order of hundreds to thou- 
sands of years. Alternatively, f tl values can be evaluated from 
evolutionary codes including shell-burning and reasonable pre- 
scriptions for mass-loss (e.g. Wood & Faulkner 1986). From 
such models, Vassiliadis & Wood (1994) find higher f t , values 
between 3 000 and 20 000 yr, but with no obvious correlation 
with mass. Moreover, the derived f tr depends critically on the 
particular phase of the pulse cycle at which the star leaves the 



AGB. Therefore, f tr is very uncertain, and even for the opti- 
mistic case in which the mass-loss rates were known in detail, 
a large star-to-star variation of f t , might be expected. 

In our models, we decided to compute PN grids for single 
values of constant f tt . This approach is useful for calibrating the 
typical values of f tr , which is otherwise essentially unknown. In 
Paper I, we found that data about Galactic PNe, as for example 
the relation between mass and radius of the ionized matter, can 
be reproduced well with shorter transition times of 1000 years 
or below. 

The type of post-AGB evolution is also determined, mainly, 
by the phase of the pulse cycle in which the star leaves the 
AGB (Iben 1984; Wood & Faulkner 1986; Vassiliadis & Wood 
1994): if it occurs during the quiescent phases between ther- 
mal pulses, a classical H-burning track is followed, otherwise 
the central star starts the post-AGB as a He-burning object. In 
principle, both H- and He-burning cases are to be expected in 
samples of PNe. The main challenge of theory, in this case, is 
to constrain the expected ratio of H- and He-burning stars. 

In our models, we compute PN grids for either H- or He- 
burning tracks. The resulting synthetic samples can then be 
mixed in varying proportions, before they are compared with 
data. 

For the present study extensive grids of PN models are cal- 
culated with the following input parameters: 

- Z = 0.004, 0.008, 0.019, which correspond rougly to the 
SMC, LMC, and solar composition; 

- ? tr = 500, and 1 000 yr; 

- both H- and He-burning tracks. 

Each grid corresponds to at least 20 tracks computed for stellar 
masses between M; = 0.6 and 5 M Q . All models adopt the f(L) 
prescription for the terminal wind velocity (Paper I). Additional 
sets with different f te values - namely 800 and 2500 yr - are 
calculated for the Z = 0.019, H-burning case. 

3. Properties of synthetic PN samples 

3.1. Modelling populations of PNe 

In order to simulate populations of PNe in external galaxies, 
we employ a simple population synthesis approach. Since the 
PN lifetimes, Af PN , are much shorter than the main sequence 
lifetimes of their progenitors, t h , the initial mass distribution 
of PNe - the number of PNe from stars with initial mass M; - 
is fairly well described by: 

N(M0 oc 0(Afi) i//(t - t h ) Af PN , (3) 

where <p(Mi) is the initial mass function (IMF), and t/f(f) is the 
star formation rate (SFR, in mass per unit time) at the galactic 
age t. 

In the present work, we simply assume a Salpeter (1955) 
IMF, 4>{M) oc Mr 2 35 . The stellar quantities, T H (Mj,Z) and 
Af PN (Mi,Z), are interpolated from the existing grid of models 
to any given value of initial mass M; and metallicity Z. Given 
an age t, the stellar metallicity is selected from a given age- 
metallicity relation (AMR) Z(t). 
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By assuming a given star formation history \]/(t) and Z(f), 
Eq. Q is used to construct the probability distribution function 
of initial masses. Then, any synthetic sample of PNe is ran- 
domly generated with the aid of a simple Monte Carlo method, 
provided that the above distribution is complied with. Hence, 
for each star, a random age is selected in the interval [0, Af PN ]. 
Finally, all stellar parameters (nebular masses, radii, velocities, 
emission fluxes, and the properties of central stars) are obtained 
from interpolation in the existing PN model grids. 

It should be recalled that the resulting synthetic samples 
strictly obey the basic theories of stellar evolution and popu- 
lation synthesis of galaxies, and therefore they present some 
of the intrinsic selection effects and biases that are inherent to 
these theories. For instance, tracks with higher PN-phase life- 
times are naturally more represented in these samples. In or- 
der to simulate observed samples, additional selection criteria 
should be imposed, such as e.g. a given limit in [Om] /15007 
magnitude and/or surface brightness. In the present work, such 
selection criteria are assumed a posteriori, when required for 
the comparison with empirical data. 

For each one of PN model grids, we simulate synthetic PN 
samples until a number of 10 4 PNe brighter than M5007 — is 
reached. This guarantees an adequate sampling of the brightest 
part of the PNLF. 

3.2. Comparison with empirical data 



Figures[3H6|illustrate the predicted properties of synthetic sam- 
ples of PNe and their central stars for exemplifying choices of 
the input parameters (namely: constant metallicity Z = 0.019; 
ftr = 500 yr; H/He-burning tracks; constant SFR) and selection 
criteria (specified in the captions). Similar plots can be obtained 
for different input parameters. 

Let us start with Fig. [5] that shows the predicted distribu- 
tion of the PN central stars in the HR diagram, as a function of 
the [O m] /I5007 luminosity (larger circles correspond to more 
luminous PNe). On average, central stars with higher lumi- 
nosity and larger effective temperature produce brighter PNe 
in the [Om]/I5007 line. We see that PNe - emitting in the 
[O m] A5Q01 line - should be optically-thin to the H-ionising 
flux as long as their central stars are evolving along the hor- 
izontal part of their post-AGB tracks; later they may become 
optically-thick as the PN nuclei start fading. Moreover, we no- 
tice that He-burning central stars of optically-thick PNe tend to 
concentrate towards lower luminosities compared to H-burning 
objects. 

Figures |4j|6] display how synthetic samples of PNe dis- 
tribute over a few relevant diagrams, where observed PNe are 
found to follow clear correlations between the involved nebular 
parameters, namely: the ionised mass, the nebular radius, the 
electron density and specific intensity line ratios. For a com- 
plete discussion of the diagrams we refer to Paper I, where the 
empirical data are compared to predictions of individual pa- 
rameters as a function of the PN age for a representative set 
of PN models. In our present study, the comparison is instead 
performed in terms of statistical properties, with the aid of syn- 
thetic populations of PNe. 
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Fig. 3. Location of PN central stars in the HR diagram. The 
simulated sample corresponds to PNe with M5007 < 1 out of 
a total population of 500 objects, assuming the following set 
of input parameters: constant SFR, constant metallicity Z = 
0.019, transition time f tr = 500 yr, H- or He-burning post-AGB 
evolutionary tracks. Filled and empty circles denote PNe that 
are optically thick and optically thin to the hydrogen ionising 
photons, respectively. Circles are larger for brighter M5007 




0.001 



Fig. 4. Ionised mass-radius relation for Galactic PNe. 
Observed data (crosses) are taken from Zhang (1995) and Boffi 
& Stanghellini (1994). The synthetic sample (circles) is the 
same one as in Fig. [3] 
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Fig. 5. Electron density-radius relation for Galactic PNe. 
Observed data (crosses) are taken from Phillips (1998). The 
synthetic sample (circles) is the same one as in Fig. [21 



Fig. 7. Anticorrelation between He n A4686 and [O ra] ,15007 
lines. The same as Fig. [6] but for Z = 0.008 



well, which confirms the good agreement already obtained in 
A few remarks should be made at this point. First, we Paper I. 
see that the synthetic samples recover the observed data quite 
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I([OIII]S007) / I(HeII4686) 



Fig. 6. Anticorrelation between Heii/14686 and [Om] ,15007 
line intensities. Observed data (crosses) are taken from the 
Strasbourg-ESO Catalogue of Galactic Planetary Nebulae 
(Acker et al. 1992). The synthetic sample (circles) is the same 
one as in Fig. [3] 



7500 10" 1.25x10" 1.5x10" 1.75x10" 2x10" 

Fig. 8. Distributions of electron temperatures in PNe. Shaded 
histogram: observed distribution of T s [0 in] for a sample of 99 
Galactic PNe taken from McKenna et al. (1996). Empty his- 
tograms: predicted distributions derived from synthetic sam- 
ples with parameters Z = 0.019, f tl = 500 yr, H-burning tracks 
(continuous line) and He-burning tracks (dot-dash line) 
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100 x I([OIII]5007)/I(Hj8) 

Fig. 9. Distributions of [O in] A5001 intensity relative to H/3, 
on the scale I(H/3) = 100. Shaded histogram: observed dis- 
tribution for 886 Galactic PNe taken from the ESO-Strasbourg 
Catalogue (Acker et al. 1992). Continuous-line histogram: pre- 
dicted distribution derived from synthetic samples with param- 
eters Z = 0.019, f t , = 500 yr, H-burning tracks, re-scaled to 
the total number of observed PNe. Dot-dashed line histogram: 
the same but for He-burning tracks. We remark that the cor- 
responding distributions for t a = 1000 yr change little with 
respect to those already plotted 

For the ionised mass-radius and electron density - radius 
relationships (Figs.|4]and|5} we have the interesting result that, 
for the given set of adopted parameters, the observed excursion 
in radius is accounted for by including both H- and He-burning 
tracks, the latter group mainly contributing at larger radii. This 
is explained by the larger evolutionary timescales of the He- 
burning stars, so that more expanded nebulae are expected at 
given ionised mass. We also note that the observed points at 
lower masses in Figs. 0] could be recovered by widening the 
selection criteria in our simulated sample, i.e. including also 
PNe that are fainter than M5007 = 1 (which is our present upper 
limit in magnitude). 

The empirical anticorrelation between the Heii/14686 and 
[Om]/I5007 line intensities, displayed in Fig. |6] is well re- 
produced by the synthetic sample for metallicity Z = 0.019, 
which however follows a narrower relationship than observed. 
A larger dispersion is obtained by a simulated PN sample with 
the same parameters but for a lower metallicity, i.e. Z = 0.008 
(see Fig.0. Intepreting this result in terms of the oxygen abun- 
dance X(0), we may derive the indication that Galactic PNe 
derive from stellar progenitors with lower X(0) than assumed 
in the solar-metallicity models (Anders & Grevesse 1989). To 
this respect, it is worth recalling that, according to the recent 
determination by Allende Prieto et al. (2001), the Sun's oxy- 
gen abundance is lower by about 0.2 dex than previous esti- 



mates (e.g. Anders & Grevesse 1989), that are actually adopted 
in our solar-metallicity models. 

Finally, Figs. [8] and |9] show the comparison between 
observed and predicted distributions of two relevant quan- 
tities, namely: the electron temperature and the line ratio 
7([Oin]/15007)//(H/3), on the scale /(H/3) = 100. The empir- 
ical data refer to samples of Galactic PNe. 

As already mentioned, the electron temperature is a criti- 
cal quantity when dealing with the intensities of collisionally- 
excited lines, such as [Oin] A50Q1. It follows that the distribu- 
tion of r e represents an important observable to test our PN 
models, particularly in relation with the [O in] /15007 PNLF. 
We see in Fig. [8] that the comparison turns out successfully 
as models are able to recover the location of the distribution 
peak, centered at around T e = 10000 K. We also notice that 
the the H-burning set of models appears to describe well all the 
main features of the observed T e histogram, the He-burning set 
populating a narrower T e range than observed. 

The distribution of the [O m] A5QQ1 intensity relative to H/3 
also represents a relevant constraint to model predictions. It 
is worth remarking that a satisfactory reproduction of the ob- 
served data is naturally obtained with our PN models - mostly 
the H-burning set - without introducing any additional assump- 
tion or free parameter. Again, this sets another important im- 
provement on previous studies of the PNLF (e.g. Mendez & 
Soffner 1997), in which the reproduction of the /([Oin] ,15007) 
distribution - without the use of a photoionisation code - is de- 
rived analytically and it requires a number of calibration steps, 
i.e. the /([Oin] /15007) values are adjusted during the simula- 
tions in a fashion depending on the mass of the central star and 
its position in the H-R diagram. 

4. Properties of synthetic PNLFs 

In this section, we will discuss the dependecy of the PNLF cut- 
off on several parameters. With the aid of our models, we start 
answering the most basic question, namely which stars produce 
the PNe populating the cut-off (Sect. 14. lV Then, we analyse the 
intrinsic dependencies on two poorly known model parameters, 
namely the dominant nuclear burning regime of post- AGB cen- 
tral stars (Sect. l4~2l and the transition time (Sect. l43t . Finally, 
we move to investigate the effect of two genuine parameters 
of stellar populations - that are likely to vary from galaxy to 
galaxy -, namely age and metallicity (or rather the oxygen 
abundance) of PN progenitors (Sects. I4~4l and l4.5l respectively). 

4.1. Which PNe form the cut-off? 

This simple and basic question has received different answers 
by different authors. Just to mention two of them, Mendez et al. 
(1993) conclude that the cut-off PNe are optically-thin objects, 
whereas Dopita et al. (1992) derive that some of the brightest 
LMC PNe are optically thick. Mendez et al. (1993), Mendez & 
Soffner (1997), and Dopita et al. (1992) conclude that the cut- 
off PNe are young objects, at least in star-forming galaxies like 
the LMC (i.e. ~ 0.8 Gyr), whereas Ciardullo & Jacoby (1999) 
argue that, due to circumstellar extinction mainly affecting the 
most massive AGB stellar progenitors, the brightest observed 
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PNe in galaxies have always higher ages (i.e. > 1 Gyr). In our 
models, we attempt to answer this question just by analysing 
the brightest PN models in our available grids, without making 
any additional assumption about their properties. 

We start by considering Fig. E3 that displays the run of 
M5007 versus age, for the grid of Z = 0.008, f tt = 500 
yr, H-burning tracks. This grid was chosen mainly because it 
presents, for the case of continuous SFR, a cut-off magnitude 
M* = -4.63 (see Fig. [H] continuous-line histogram), very 
close to the observationally-calibrated value of Ciardullo et al. 
(2002). From Fig.^|it turns out that this cut-off M* coincides 
with the maximum M5007 reached by the 2.5 M track, which 
defines the characteristic progenitor stellar mass of cut-off PNe. 
Similar results are obtained for each grid of tracks; what varies 
from case to case is just this characteristic mass, which tends to 
be slightly lower for higher thansition times and for He-burning 
tracks. Anyway, in all cases this characteristic mass is com- 
prised between 2 and 3 M Q . 

What determines the [O in] A50Q1 maximum luminosity of 
the individual tracks ? For most of them, this maximum occurs 
when the PN switches from optically thin to optically thick, 
i.e. when yUe changes from < 1 to = 1. In fact, at that point 
the central star starts dimming towards its WD cooling track, 
no longer providing enough photons to completely ionize the 
nebula, so that the Stromgren sphere suddenly starts decreas- 
ing. At the same time the ++ region also starts receeding in 
radius, and hence to the [Oin] A5Q07 luminosity decreases. 

From these considerations it follows that the PNe observed 
close to their maximum [O ra] A50Q1 luminosity, can be either 
optically thin or thick, depending on whether they are evolv- 
ing along the ascending or descending part of their tracks in 
Fig-Ell This implies that both thin and thick nebulae should be 
present in the PNLF cut-off, the frequency of optically-thin ob- 
jects likely being higher given the faster evolution evolution of 
M5007 in the descending branch. In this detail, our results differ 
from the assumptions/conclusions of both Mendez et al. (1993) 
and Dopita et al. (1992). 

At this point it is natural to wonder: Why is the absolute 
M5007 minimum value attained just by the 2.5 M Q tracks ? The 
answer to this question is not trivial, since a lot of factors are 
envolved. Our proposed explanation is based on the temporal 
behaviour of three quantities, that are plotted in Fig.^2 as a 
function of the PN age for a few models with CSPN masses 
around 2.5 M Q . They are: 

I) The ionisation rate for ++ , that expresses the number of 
available ionising photons emitted by the central star per unit 
time. This rate is given by 4nR 2 Q(0 ++ ), where R is the stellar 
radius and Q(0 ++ ) is the ionising flux. We recall that, due to 
the similarity between their ionisation potentials, ++ and He + 
have similar ionisation fluxes Q (see Paper I). Therefore the 
behaviour of the ionisation rate for ++ mimics that of He + , 
which is displayed in Fig. ^2 We see that larger ionisation rates 
correspond to more massive stellar progenitors, as they produce 
more massive and luminous central stars. Therefore, models 
with progenitor masses of 3.0 and 3.5 M Q could be potentially 
brighter in [Om] /15007 than the 2.5 M Q model. 

II) The maximum nebular mass available to the ++ ionis- 
ing photons. This is given, at each instant, by the mass of the 
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Fig. 12. Theoretical IFMR adopted in the present calculations. 
Notice that at given initial masses at decreasing metallicity 
models predict higher final masses, hence more massive cen- 
tral stars of PNe. 

shocked shell, M s h, which is steadily increasing with time. It 
is clear that, in principle, a higher [O in] /15007 luminosity is 
favoured at larger M s ^. 

Ill) The electron density, iV e . In fact, it is well-known that 
electrons do play the major role in populating by collisions the 
upper level of [O m] /15007 so that the higher N s the stronger 
the spectral line. This holds until a critical density A^ nt is 
reached, above which the collisional de-excitation becomes im- 
portant and the the intensity of [O m] /15007 starts to weaken. 
Moreover, we see that N e steadily decreases with time. 

The answer to the initial question rests on the behaviour 
of all above quantities, whose evolution is controlled by the 
timescales of both central star and nebular dynamics. From 
Fig-^Dit turns out that during the initial PN evolution, i.e. in- 
cluding the stages that bracket the ionisation rate's maximum, 
only a fraction of the hydrogen ionising photons are effectively 
absorbed by the nebula that will be optically thin (//h < 1), 
while a large fraction will contain ++ . The comparison be- 
tween the panels for different masses in Fig. ^]then reveals 
the following. 

On one hand, for the most massive tracks (k, 3 M Q ) the 
evolutionary rate of the central star is so rapid that, during the 
stages of significant ionisation rate, the shocked shell has en- 
gulfed a relatively little amount of mass. It means that, despite 
the large availability of energetic photons to ionise ++ , this 
is effectively reduced by the little availability of nebular mat- 
ter. Additionally, we notice that at these ages a further weak- 
ening of [O in] /15007 emission may be produced as the neb- 
ular densities approach, or even overcome, the critical limit 
for collisional de-excitation of the [O in] /15007 upper level, 
i.e. ~ 4.9 x 10 5 cirT 3 . All these factors concur to prevent PN 
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Z = 0.008, SFR = const., t tr =500 yr, H-burn. 




PN age (10 3 yr) 

Fig. 10. Plot of M5007 versus PN age for the tracks in our grid of Z = 0.008, t a = 500 yr, H-burning tracks. It can be noticed that 
the most luminous PNe are those of mass close to 2.5 M . Optically thick nebulae are marked with continuous lines, optically 
thin with dotted ones. 
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Fig. 11. Several model 
quantities plotted as a 
function of PN age, for 
initial masses between 2 
and 3.5 M Q , in our grid of 
Z = 0.008, ttr = 500 yr, 
H-burning tracks. The age 
interval covers the com- 
plete period in which the 
object is seen as a bright 
PN in M 50 07. Optically- 
thick nebulae are marked 
with continuous line, those 
optically-thin with dotted 
line. The 2.5-M© track, 
responsible for the PNLF 
cut-off, is the most massive 
one for which the periods 
of high photoionisation 
rate 47rfl 2 g(He + ) partially 
coincide with a significant 
shell mass M s h, and with 
an electron density 7V e 
below the critical value for 
collisional de-excitation of 
the [Om] ,15007 line. For 
lower-mass stars, the shell 
masses become higher, but 
the photoionisation rates 
much lower 
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models, having progenitor stars with M > 2.5M , from be- 
coming brighter in [O m] A5QQ1 than those produced by stars 
with M ~ 2.5M . 

On the other hand, for less massive tracks (< 2.5 M ) there 
is always a coincidence of both significant ionisation rates and 
relatively large nebular masses, due to the slower evolutionary 
rate of the central star. However, the corresponding PNe are 
still fainter in [Oin] /15007 than the 2.5 M model because of 
the lower ionisation rates. 

In summary, PN tracks with progenitors masses of ~ 2.5 
M present the right combination of the dynamical timescale 
of the nebular expansion (that determines M s h and N e ), and 
the evolutionary timescale of the central star (that determines 
47i7? 2 g(He + )), to maximize the [O ra] ,15007 luminosity. 

An additional effect in determining the PNLF cut-off is re- 
lated to the PN lifetime. PNe produced by more massive pro- 
genitors stars are characterised by shorter Af PN and therefore, 
according to Eq. Q, they contribute with fewer objects to the 
PNLF. In this regard, it should be noticed that, in their very 
early M5007 evolution, some of the PN tracks of higher ini- 
tial masses (k. 4 M ) may present spikes as bright as M5007 = 
-6.5 2 . However, these spikes do not affect at all the position of 
the PNLF cut-off because they are too short-lived (^ 100 yr). 

Finally, we recall that 2.5 M tracks correspond to main se- 
quence lifetimes of just 7 x 10 8 yr. The minimum M5007 magni- 
tudes become progressively fainter as we go to tracks of lower 
masses (larger turn-off ages). These facts have immediate im- 
plications for the dependence of M* on the galaxy's star forma- 
tion history (Sect. 03). 

The next related question is now: What are the masses of 
the CSPNe forming the PNLF cut-off? According to our AGB 
calculations (Marigo et al. 1999) a star with initial mass of 
~ 2.5 M Q and metallicity Z = 0.008 ends its evolution as 
bare carbon-oxygen core of about ~ 0.75 M Q (see Fig. I12> . 
In Fig.[0]we show how the mass distributions of the CSPNe 
are predicted to vary if a limit minimum [O in] /15007 magni- 
tude M is applied to the synthetic samples, i.e. we consider 
only PNe brighter than M. We see that at decreasing M, the 
peak of the theoretical distribution shifts towards larger masses, 
for less massive and fainter CSPNe fall outside the adopted lu- 
minosity cut. As expected, in the extreme case of M — -4.5 
(close to the observed PNLF cut-off) the peak is centered just 
at Mcspn ~ 0.75 M (corresponding to the 2.5 M progeni- 
tor star). An interesting point to be noticed is that, contrarily 
to what often assumed, the CSPN mass distribution appears to 
differ from that of white dwafs (bottom-right panel of Fig.ll2i. 
The latter is in fact characterised by a sizeable tail towards 
larger masses, which is not expected in the case of CSPNe. 
This can be understood by considering that the evolutionary 
speed of white dwarfs slows down so much that their lifetimes 
can be considered practically infinite (regardless the mass of 
the stellar progenitors), which determines a sort of accumula- 
tion of objects in the high-mass bins. Other differences between 
the mass distributions of CSPN and white dwarfs may be intro- 
duced by observational biases, e.g. selection effects against the 

2 These spikes show up mainly in the sets of PN models with f tt = 
1000 yr. 



Z = 0.008, SFR=const., t lr =500 yr 



_H-burn., all ages 
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M(5007) 

Fig. 15. Effect of the nuclear burning regime of central stars 
(H- or He-burners) on the PNLF. The histograms show PNLFs 
as derived from simply counting PNe in 0.2-mag bins, whereas 
the smooth lines show the result of fitting Eq. (|2j to the bright- 
est 1-mag interval. For PN samples containing stars of all 
ages, the cut-off is affected just slightly by the nuclear burning 
regime of central stars - in fact in the case here illustrated M* 
is -4.63 for the H-burning case, and -4.60 for the He-burning 
one. Only when we limit to consider stars of older ages, the na- 
ture of the tracks is of more relevance, He-burning PNe becom- 
ing much fainter than H-burning ones: after 2 Gyr, M* becomes 
-3.64 and -1.96 for H- and He-burning tracks, respectively. 



detection of PNe with slowly evolving, low-mass, central stars 
(Stasihska et al. 1997; see also Phillips 2000 for a critical anal- 
ysis of the empirical method used to evaluate N(Mcspn))- 

4.2. Dependence on H/He burning tracks 

As discussed in Sect. 12.41 the relative proportions of PN H- 
or He-burning nuclei is essentially unknown. An estimation is 
given by Dopita et al. (1997) who derive a frequency ratio He/H 
~ 45% - 70%, basing on evolutionary timescales and nebular 
sizes for a sample of 15 LMC PNe. 

In general, He-burning tracks are characterized by longer 
evolutionary timescales and fainter luminosities than those of 
H-burning models, (see also Paper I). The resulting evolution 
of M5007 is illustrated in Fig. \^\ that should be compared to 
the corresponding figure for H-burning tracks of Fig. EH It is 
evident that He-burning PNe evolve slowlier, and preferentially 
in optically thin conditions. However, the characteristic masses 
of the PNe at the cut-off change only slightly. 

Figure [TBI shows a direct comparison between PNLFs de- 
rived from H- and He-burning tracks, assuming Z = 0.008, 
f tr = 500 yr, and continuous SFR. The shape of the PNLF 
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Fig. 13. Theoretical mass 
distributions for initial 
metallicity Z = 0.008. 
The left and top-right 
panels refer to central stars 
of PNe more luminous 
than a limit magnitude 
M5007 as indicated, with 
continuous and dashed 
lines corresponding to H- 
and He-burning stars, re- 
spectively. For the sake of 
comparison, the predicted 
mass distribution of white 
dwarfs is shown in the 
bottom-right panel 
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Fig. 14. The same as Fig.^| but for He-burning tracks. Notice also the wider range of PN ages compared to the H-burning set 
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Fig. 16. Effect of the transition time on the PNLF, exemplified 
with the aid of the Z = 0.008, H-burning models. For PN sam- 
ples containing stars of all ages, the cut-off is affected signif- 
icantly by f t , - in fact M* changes from -4.63 to -4.33 as f tt 
varies from 500 to 1000 yr, in the first two cases here illus- 
trated. When we limit only to stars of older ages, the effect of 
t u becomes of little importance. 



somewhat changes, but not that dramatically, so that a fitting 
of Eq. (0 to the brightest PNLF bins is straightforward. In fact, 
an unweighted least-squares fitting of Eq. [5]- performed on a 
PNLF with 0.2-mag wide bins, and limited to the brightest 1- 
mag interval - produces M* = -4.63 for the H-burning case, 
and M* = -4.60 for the He-burning one, i.e. yielding almost 
identical values. Similar results are obtained by using different 
transition times and metallicities, while adopting either H- or 
He-burning tracks. 

However, we note the M* tends to decrease much faster 
with age in the case of He-burning tracks, than H-burning ones. 
This effect can be appreciated in Fig.^] and is discussed fur- 
ther in Sect.l4"4"lbelow. 

4.3. Dependence on the transition time 

Whereas there is no reason to believe that this is the case, the 
lack of robust predictions for the transition time (see Sect. 12. 4> 
as well as practical reasons impose that we first simulate 
PNLFs under the assumption f tr is the same for all stars. In 
fact, when we do this with our available tracks (for f tl - = 500 
and 1000 yr), we obtain the results depicted in Fig.^] a high 
dependence of the cut-off magnitude M* on the transition time. 
This can be explained by the large impact of f tr on the timescale 
for nebular expansion, especially for PNe of higher masses, 
whose evolutionary timescales are comparable to f tr itself. 



Fig. 17. The same as Fig.[l6] but for Z = 0.019 and 4 different 
values of f tl , and limited to the case of continuous SFR. The 
trend of M* with f tl - is very regular, with values M* = -4.81, 
-4.74, -4.67, and -4.23 for = 500, 800, 1000 and 2500 yr, 
respectively. 



In order to investigate further this behaviour, we computed 
additional sets of models for Z = 0.019, H-burning tracks, and 
f tr = 800 and 2500 yr (Fig. I17> . If we consider PN samples 
that cover all ages, the cut-off exhibits a very regular trend, 
becoming brighter at shorter f^. The maximum [Om]/i5007 
brightness a PN sample may attain, however, will be limited 
by the fact that f tr cannot be shorter than a few hundreds of 
years. The value of f tl = 500 yr, succesfully tested in Paper I 
and in this work on a number of Galactic PN observables (see 
Sect. 13. 21 . is already an extremely low value if compared to 
present results of evolutionary calculations (e.g. Vassiliadis & 
Wood 1994). Actually, it yields M* values that are slightly too 
bright, i.e. ~ -4.8 for f tr = 500 yr, than the observationally- 
calibrated M* = -4.47, while somewhat lower M* values are 
obtained with f tl = 800 yr and f tl = 1000 yr. On the other hand, 
f tr longer than 2500 yr would produce too faint M* , i.e. -4.3. 
In summary, according to our present PN models, the observed 
cut-off magnitudes M*, for galaxies with ongoing SFR, will be 
obtained for f tr values between 500 and 1000 yr. 

At this point two important remarks are worthy to be made. 
First, as far as the bright PNLF cut-off is concerned, the indi- 
cations we derive on f tr should be taken as a sort of calibration 
that applies only to those PNe that actually form the cut-off, 
namely PNe with relatively massive central nuclei (Mcspn ~ 
0.7 - 0.8M o ), produced by relatively massive (M ~ 2.5M ) 
and young (t h ~ 0.7 Gyr) stellar progenitors. In this range our 
results agree with the predictions for f tr derived by Stanghellini 
& Renzini (2000) from analytic relations (see their figure 11). 
It is clear that, more realistically, t a > 500 - 1 000 yr should be 
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Fig. 18. Effect of the age distribution on the PNLF. To do this 
plot, we have generated a synthetic sample of PNe from a con- 
stant SFR, and then considered in the PNLF only the PNe with 
ages older than the specified age limit. It may be seen that the 
cut-off depends very much on the age of the latest star forma- 
tion episode, and hence on the recent (less than a few Gyr) star 
formation history. 



t tr =1000 yr, H-burn. . 

t lr =1000 yr, He-burn. . 

t tr =500 yr, H-burn. 
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Fig. 19. Variation of the brightest magnitude bin with the age 
of last SFR. This bin behaves similarly to M* (although al- 
ways 0.2 - 0.4 mag fainter), but is more easily defined for all 
ages. The curve has been derived by searching for the steep- 
est change in the logarithmic of number counts, in a series 
of PNLFs generated in the same way as in Fig. [D^] Different 
cases (H-/He-burners, two transition times) are shown, always 
forZ = 0.008. 



associated to PNe with lower-mass stellar progenitors, but this 
would not affect at all our results on the cut-off brightness. 

Second, a further support to our estimated ranges for f tr and 
PN phase duration (see for instance Figs. 1101 and [T4T > comes 
from the kinematic ages (from less of 1 000 yr up to almost 
5 000 yr) empirically derived by Dopita et al. (1996) for a HST 
sample of 15 PNe in the LMC. Many of these PNe are, in fact, 
classified as Peimbert type I PNe, suggesting that they likely 
hide massive central stars. 

4.4. Dependence on age of the last burst of star 
formation 

The age of PN progenitors has long been suspected to have 
a major effect on the cut-off magnitude (Dopita et al. 1992; 
Mendez et al. 1993; Stanghellini 1995; Jacoby 1996). Age is 
closely related to the turn-off mass, and to the remannt mass 
via the IFMR. For any reasonable relation of this kind, at older 
ages only PN central stars of lower mass, and hence of lower 
intrinsic luminosity, should form. Our computations fully con- 
firm the expected trends. 

To illustrate the effect of the age, we simulate PN samples 
adopting a SFR which is constant over the age range [?i ast , f ga i]. 
The upper limit, f ga i, denotes the age of the oldest stars in a 
galaxy, and here it is assumed constant and equal to 12 Gyr. 
The lower limit, fi ast , is the age of the most recent burst of star 
formation, and it is let vary between - for systems with on- 



going star formation, namely the disk of spirals and irregular 
galaxies - and f ga i - for predominantly old systems, like haloes 
and bulges, and old ellipticals. 

The results for the PNLFs are illustrated in Fig. We 
clearly see that, as /i ast increases, the PNLF becomes progres- 
sively depleted of its brightest PNe, and the cut-off shifts to 
fainter magnitudes. At ages older than a few Gyr, the PNLF 
may also change its the characteristic morphology, i.e. shaping 
like a plateau with a sharp cut-off in place of the smoother trend 
that is seen for younger ?i ast . 

Figure[23shows the behaviour of the PNLF over the bright- 
est magnitude bin as a function of age. This bin is simply de- 
fined as the one of the steepest change in the logarithm of num- 
ber counts. It behaves similarly to the cut-off M*, but it is sys- 
tematically fainter than that quantity by 0.2 to 0.4 mag 3 . As a 
matter of fact, the fading of the brightest magnitude bin starts 
just after fi ast becomes older than the age of the brightest PN 
progenitors - that means after ~ 0.7 Gyr, which is the turn-off 
age of ~ 2.5 M stars. Then, at older fi ast , the dimming of the 
cut-off is clearly seen for all ages. 

Finally, we remark that on the base of our H-burning PN 
models, the fading rate of M* with fi ast is of roughly 0.5 

3 The reason why we plot the brightest magnitude bin, is that it 
behaves in a much smoother way than M* , for all values of fi ast . In 
fact, for f las , older than say 4 Gyr, the PNLF often presents too sharp a 
cut-off followed by a shallow plateau, and a good fitting of Eq. is 
not possible. 
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Fig. 20. Effect of mean metallicity on the PNLF. For constinu- 
ous SFR, the cut-off depends on metallicity in a non-monotonic 
way. For older populations, instead, the behaviour is mono- 
tonic. See text for details. 



mag/Gyr. This rate is higher than those predicted by other au- 
thors (e.g. Jacoby 1996) with the aid of simpler models. 

4.5. Dependence on metallicity 

Figure |20| illustrates how the PNLF is predicted to vary with 
metallicity according to our present PN models, for the cases 
of fiast = and fi ast = 2 Gyr. It can be noticed that our models 
exhibit an appreciable dependence of M* on metallicity. In the 
case of continued SFR, this dependence is not even monotonic 
with metallicity, being the Z = 0.008 cut-off ~ 0.4 mag fainter 
than those found for Z = 0.004 and 0.019. For f last > 2 Gyr, 
instead, we find that more metal-rich PN samples correspond 
to fainter cut-offs. The same behaviour is shared by all our sets 
of models. 

These results may seem somewhat unexpected, consider- 
ing the previous claims for a weak dependence of the cut-off 
on metallicity (e.g. Ciardullo et al. 2002). To understand this 
point, first it is worth briefly recalling the origin of these claims. 
The most comprehensive analysis of the [O/H] dependence is 
carried out by Dopita et al. (1992) on the base of photoionisa- 
tion models of optically-thick nebulae. It is found that the peak 
conversion efficiency of stellar photons to nebular [O in] /15007 
ones is given by 



l °s\ F },l-l \ = -0.8791 +0.1459 [O/H] 
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Fig. 21. Variation of M* with [O/H], according to the quadratic 
relation (Eq. |5j continuous line) derived by Ciardullo et al. 
(2002) on the base of Dopita et al.'s (1992) data (circles). Filled 
circles rely on the assumption that the stellar progenitors of the 
cut-pff PNe are 0.8 Gyr old, whereas empty circles show the 
location of the LMC and M31 if their SFR stopped 5 Gyr ago. 
The adopted [O/H] values are taken from Dopita et al. (1992) 

for the absolute luminosity of the star at the PNLF cut-off in 
a few galaxy systems (table 2 in Dopita et al. 1992), one can 
derive the following relation, quoted by Ciardullo et al. (2002): 



5M* = 0.928 [O/H] + 0.225 [O/H] + 0.014 



(5) 



where [O/H] is the oxygen content in the nebula scaled to the 
Sun's. By using this formula, together with the inferred values 



which has indeed a weak dependence on the oxygen abun- 
dance, as displayed in Fig.|^(filled circles). We should how- 
ever emphasize that this latter relationship holds only under a 
precise condition, namely the LMC, SMC, and Galactic Bulge 
contain young and massive PN progenitors, of ages similar to 
0.8 Gyr. If young stars were not present in the galaxies stud- 
ied by Dopita et al. (1992), the fitting relation by Ciardullo et 
al. (2002) would not be valid, as it can easily be understood by 
looking at the empty circles in the figure. Another interesting 
point is that, according to Eq. (Q, M* should reach a minimum 
for near-solar metallicities, contrarily to our findings. 

Following Ciardullo (2003), the results of Fig. |2] are the 
consequence of two competing factors: (1) since oxygen is one 
of the main coolants of a nebula, the flux in [O in] /15007 in- 
creases only with the square root of the oxygen content (Jacoby 
1989); and (2) the UV flux coming from the CSPN decreases 
with metallicity. This latter trend, actually, is dictated by the 
earlier onset of the superwind, and hence the lower final stellar 
masses and luminosities, that are expected in AGB models of 
higher metallicities. 

After these considerations, how can we understand our re- 
sults of Fig. [20] in terms of basic stellar parameters? First of 
all, in agreement with other studies (see e.g. Lattanzio 1986; 
Vassiliadis & Wood 1993; Wagenhuber & Groenewegen 1998; 
Marigo 2001), our underlying AGB models (Marigo et al. 
1999) present a clear trend of decreasing final masses at higher 
metallicities, as depicted in Fig. El This behaviour ensures 
that, for given stellar age and initial mass, lower-metallicity 
PNe will have central stars with higher masses and hence UV 
fluxes, in agreement with Ciardullo's reasoning. We also re- 
mark that although the shape of the IFMR depends on a num- 
ber of model details - such as the presence of convective core 
overshooting in previous evolutionary phases, the number and 
efficiency of third drege-up events, etc. - the key-role should 
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Fig. 22. Oxygen abundance at the end of the TP-AGB evolution 
for the tracks used in this work, as a function of the core mass. 



be ascribed to the adopted prescription for mass loss as shown, 
for instance, by Groenewegen & de Jong (1994). 

The second important effect to be considered is the change 
of oxygen content [O/H] with iron content [Fe/H]. The ini- 
tial main-sequence [O/H] values in our stellar models are sim- 
ply scaled to the Sun's determination 4 according to Anders 
& Grevesse (1989), thus yielding [O/H] = (-0.7,-0.4,0) for 
Z = (0.004, 0.008, 0.019), respectively. 

It is worth recalling that stellar evolution models predict the 
surface oxygen abundance remains essentially unaltered after 
the first and second dredge-up episodes (see e.g. Girardi et al. 
2000), while it may be affected to some extent during the TP- 
AGB phase by the third dredge-up, and possibly by hot-bottom 
burning. 

This latter process takes place in the most luminous and 
massive TP-AGB stars (with M 4M ) and it may deplete 
the stellar surface of oxygen (via the CNO cycle) in the case 
of high efficiency and sufficiently long duration of the CNO 
nuclear reactions. This circumstance is not met in any of the 
TP-AGB models here considered. 

On the contrary, the third dredge-up may cause the surface 
oxygen content to increase with time, together with the abun- 
dances of other elements like carbon and helium. The degree of 
oxygen enrichemnt in the envelope basically depends on three 
factors, namely: the efficiency of the dredge-up, the chemical 
composition of the dredged-up material, and the total number 
of dredge-up events experienced during the TP-AGB phase. 
Referring the reader to Marigo et al. (1999) and Marigo (2001) 
for more details about the TP-AGB models employed in this 
work, we recall here that 1) the efficiency of the third dredge- 
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Fig. 23. Expected variation of M* with the initial main- 
sequence [O/H] values, according to our sets of models (with 
parameters f tl = 500 and 1000 yr, and H-/He-burning CSPNe), 
as indicated (empty squares). We recall that the [O/H] ratio 
predicted in our PN models may be quite different from the 
assumed main-sequence values. See Fig. 1221 and text for more 
details 

up is parameterised by a quantity A, calibrated to reproduce 
the luminosity functions of carbon stars in the LMC and SMC; 
the oxygen abundance in the dredged-up material is assumed 
of just X{ 16 0) ~ 0.02 (in mass fraction), that represents the 
standard order of magnitude obtained by most nucleosynthe- 
sis calculations of thermal pulses (e.g. Boothroyd & Sackmann 
1988; Karakas et al. 2002; but see also Herwig et al. 1997 for 
different results). 

Despite this rather low value, it turns out (see Fig. l22l > that 
the predicted [O/H] abundances at the end of the TP-AGB 
- and during most of the PN phase - may significantly dif- 
fer from the initial main-sequence values, the deviation being 
more pronounced in models of lower metallicities because of 
the higher dredge-up efficiency and larger number of thermal 
pulses. Moreover, we notice that [O/H] does not vary at all at 
lower initial masses, say M % 1.2-1.5 M (depending on metal- 
licity), since the third dredge-up is not predicted to take place 
in those models. 

Figure l22l clearly shows that the assumption of [O/H] = 
[Fe/H], or that [O/H] is the same as in the Hn regions of a 
galaxy, may be not suitable for PNe - unless oxygen dredge- 
up does not occur in AGB stars, which is unlikely considering 
our present understanding of stellar evolution. Moreover, the 
change of [O/H] is not necessarily monotonic with [Fe/H]; in 
fact most of our models with Z = 0.004 present [O/H] val- 
ues higher than those found at Z = 0.008, but still lower than 
those for Z = 0.019. A monotonic behaviour of [O/H] with 
Z (and equivalently with [Fe/H]) is recovered only by models 
with lower masses - hence older progenitors' ages -, that do 
not experience the third dredge-up. 

In summary, if one accounts for the changes of [O/H] at the 
surface of stars with different masses and metallicities in the 
framework of the present stellar evolution models, the result- 
ing dependence of [O/H] (in PNe) with metallicity cannot be 
expressed by a simple scaling relation, rather it is characterised 
by a more complex behaviour as, for instance, illustrated in 

Fig. EH 

Such behaviour gives origin to the non-monotonic trend of 
M* with [O/H] displayed in Fig. [23] (empty squares) which, 
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Fig. 24. Peak M5007 magnitude for our models, for the several 
cases explored and as a function of CSPN mass. Notice the 
similarity with the trend of [O/H] with metallicity, shown in 
the previous Fig.l22l 



in turn, reflects in the non-mono tonic trend of M* with age. 
In order to illustrate it better, Fig. [23] shows the peak M5007 
magnitude for all our models, as a function of CSPN mass. The 
similarity of the curves with the behaviour of [O/H] shown in 
Fig.|22]is evident. 

To further investigate the sensitivity of M* to the actual O/H 
in PNe, we have performed a test on our set of Z = 0.019, ftr = 
500 yr, H-burning models, simulating a lower main-sequence 
oxygen abundance. In practice, we have re-calculated the pho- 
toionisation structure of these PN models providing CLOUDY 
with an input [O/H] lowered by a constant quantity of -0.24 
dex 5 , so that the [O/H] curve for Z = 0.019 in Fig. 1271 would 
be simply shifted downwards by this amount while keeping 
its shape unaltered. It follows that the new [O/H] curve for 
Z = 0.019 would come closer to that for the Z = 0.008 set 
at the typical CSPN masses (~ 0.70 - 0.75M o ) of the cut-off 
PNe (in the case of continous SFR, or equivalently for ?i ast = 0). 

With the aid of PNLF simulations for the new oxygen abun- 
dances we obtain that such decrease in [O/H] causes a dimming 
of the PNLF cut-off amounting to 0.48 mag for fi ast = 0, and to 
0.43 mag for fi ast = 2 Gyr. These numbers correspond to an in- 
cremental ratio AM*/A[0/H] of » 2 mag dex -1 , which predicts 
the right order of magnitude required to explain our results of 
Figs. |20] |^ and |24] i.e. in particular the difference in AM* 
between the PNLFs calculated at different metallicities. 



5 The choice of this value is motivated by the recent re- 
determination of the Sun's oxygen abundance by Allende Prieto et 
al. (2001) that is lower by just -0.24 dex than previously estimated by 
Anders & Grevesse (1989). 



This significant sensitivity of M5007 to [O/H] displayed by 
our models raises a few considerations of general validity. In 
fact, it appears clear that the predicted values of M5007 should 
crucially depend on: 

- The assumed value for the initial main-sequence oxygen 
content, hence the adopted solar oxygen abundance (if 
solar-scaled compositions are assumed, as usually is the 
case). We remark that the determination of the Sun's oxy- 
gen abundance has changed dramatically over the years - 
i.e. from log(0/H) = -3.07 in Anders & Grevesse (1989) 
to log(0/H) = -3.31 in AUende Prieto et al. (2001). This 
alone would cause a change of 0.24 dex in our scale of 
oxygen abundances. As mentioned above, a systematic de- 
crease of -0.24 dex in our O abundances - necessary if we 
were to bring our models to the solar-composition scale of 
Allende Prieto et al. (2001) - would cause a dimming of 
~ +0.4 mag in our M* values. 

- (2) The choice of the input TP-AGB models in the PNLF 
simulations. In particular, the predictions of [O/H] in PNe 
are importantly affected by several (often poorly known 
or improperly assumed) model details such as: inter-shell 
nucleosynthesis during thermal pulses, efficiency of the 
mixing episodes, mass-loss prescrition, assumed molecular 
opacities, etc. 

Finally, we recall that there is an additional reason that 
likely concurs to produce the differences in the predicted metal- 
licity dependence of the PNLF cut-off observed between ours 
and the one derived from Ciardullo et al. (2002) on the base 
of Dopita et al.'s (1992) work. In fact, while we find that both 
optically-thin and thick PNe should populate the PNLF cut-off, 
the analysis by Dopita et al. considers only the behaviour of 
optically thick objects. 

5. Discussion 

5.1. Summary of main dependences 

In this paper we have carried out a theoretical investigation of 
the [O m] /15007 PNLF and its sensitivity to various stellar evo- 
lution and population parameters. To this aim a detailed pho- 
toionisation code (CLOUDY) has been incorporated into our 
synthetic PN model (Paper I), so that we can rely on a sim- 
ple but quite complete picture of PN evolution, which allows a 
satisfactory reproduction of many PN observables like well- 
known correlations involving nebular radius, ionised mass, 
density, intensity line ratios, etc. 

In particular, major attention has been paid to the modelling 
of the bright cut-off of the PNLF, that is claimed to be a power- 
ful distance indicator in virtue of its observed invariance in host 
galaxy systems of different ages and metallicities. The results 
of our analysis on the PNLF cut-off can be summarised as fol- 
lows. In order of decreasing importance, the main parameters 
affecting the cut-off magnitude M* are: 

The age of the last burst of star formation, fi ast . It turns 
out that, over the entire range of old and intermediate ages 
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(say from 13 to 0.1 Gyr; see Fig. II (St . the possible bright- 
est PNLF cut-off is obtained and remains invariant as long as 
?i ast ^ 0.7 Gyr, that is if stars with initial masses <: 2.5M Q (and 
final masses ^ 0.75M Q ) are presently evolving through the PN 
phase. At older ages, for fi ast S 1 Gyr, we predict a sharp dim- 
ming of the PNLF cut-off, as it becomes populated by much 
fainter PNe evolved from lower mass progenitors. 

The transition time, f tr . This is one of the most ill-known 
quantities on the theoretical ground, mainly reflecting major 
uncertainties about stellar mass loss during the initial post- 
AGB stages. In addition, f tl should be affected by large vari- 
ations depending on stellar mass and phase of the last ther- 
mal pulse on the AGB. Given these premises, our assumption 
of constant transition time may appear too simplicistic, never- 
theless we are still able to derive a few relevant indications, 
namely: i) In the case of continuous SFR, the PNLF cut-off de- 
pends only on f t , of those CSPNe effectively populating it (with 
masses around 0.7 - 0.75 M Q ). A satisfactory reproduction of 
the observed cut-off brightness is obtained with f tl - ~ 500-1000 
yr. Conversely, the details of the PNLF cut-off are insensitive 
to f tr of fainter, slowly evolving PNe originated from lower- 
mass progenitors, ii) In the case of ^ ast <: 1 Gyr the dependence 
of the cut-off on f tr is weak. In fact, this time the cut-off PNe 
correspond to lower-mass central stars, for which the condition 
Af PN » f tr is usually met for reasonable values of the transi- 
tion time. In other words, the effect of t a - on the PNLF cut-off 
is practically smeared out by the long evolutionary timescales 
of the PN phase. 

The nuclear burning regime of the central star. 
Unfortunately, the relative contributions of H- or He-burning 
central stars to a PN population are still poorly determined. 
With the aid of PNLF simulations obtained assuming that 
all CSPNe are either H- or He-burning central stars, we find 
that variations in the the cut-off brightness remain moderate 
at younger ages, instead becoming huge at older ages. In 
general, we find that the PNLF cut-off for H-burning central 
stars is brighter than that for the He-burners, since these 
latter are characterised by slower evolutionary rates at fainter 
luminosities. Anyway, it is reasonable to think that both H- 
and He-burning central stars will be present in the observed 
galaxies, so that the brightest ones - i.e. the H-burning - may 
actually determine the observed PNLF cut-offs, as well as the 
age-sensitivity of the sample. 

The metallicity. Instead of metallicity, we should refer, more 
correctly, to the dependence on the oxygen content [O/H] . 
Oxygen actually plays a major role in the nebular cooling, and 
its abundance directly affects the intensity of the [O in] /15007 
line. The oxygen content in PNe traces the metallicity of the 
parent galaxy only assuming that the PN stellar progenitors 
keep their initial surface oxygen abundance unaltered during 
the entire evolution up to the PN ejection. In the framework of 
present stellar evolution models, some increase of oxygen may 
be produced by the third dredge-up occuring during the TP- 



AGB phase, the degree of enrichment being more pronounced 
at lower metallicities. Moreover, the predicted [O/H] in PNe 
depends on the assumed initial composition of the gas, that is 
usually assumed to be solar-scaled ([O/H] = [Fe/H]). All these 
factors yields a metallicity dependence of the cut-off that is not 
monotonic and also changes with age. 

5.2. On the claimed age invariance of the PNLF cut-off 

After this brief recapitulation of the main dependences in the 
models, let us now discuss to more detail the claimed age- 
invariance of the cut-off. For the cases with continued SFR, 
our models predict M* values that are very close to the ob- 
served one for the LMC and other similar galaxies in the Local 
Group. This, and the strong decay of M* with age, would sug- 
gest that we have a natural explanation for the cut-off position 
in star-forming galaxies, it being determined by stars with ini- 
tial mass of ~ 2.5 M G . This is also the conclusion reached by 
Dopita et al. (1992), based on a detailed comparison between 
PN models and observations in the Magellanic Clouds. 

However, considering that sizeable changes of M* are ex- 
pected depending on various model details, and mostly on the 
age of PN progenitors, how can we reconcile such predicted 
variations with the observational indications of a fairly con- 
stant PNLF cut-off brightness, common to galaxy systems of 
all types ? In fact, an important point is that elliptical galaxies 
in groups appear to share similar cut-offs with spiral galaxies. 

At this stage, it is worth recalling and discussing the dif- 
ferent interpretations to this issue that have been proposed in 
the recent literature, referring mostly to Jacoby (1996) and 
Ciardullo (2003), who provide two extended reviews on the 
PNLF as a distance indicator. Actually, the two reviews con- 
tain two different, perhaps complementary, explanations to the 
age-invariance of the cut-off brightness. 

Following Jacoby (1996, and references therein) the con- 
stancy of the PNLF cut-off, over the age range of 3 - 10 Gyr, 
should be essentially ascribed to the extremely narrow distri- 
bution of central star masses, say 0.58 ± 0.02 M . This would 
correspond to initial masses in the range 1-2 M accord- 
ing to the semi-empirical IFMR for the solar neighbourhood 
(Weidemann 1987; see also Herwig 1997). Our AGB mod- 
els for solar metallicity (see Fig. II 21 recover well this part 
of the observed IFMR, nevertheless the present PNLF simu- 
lations predict, on the contrary, large variations of the cut-off 
over the above age range. Actually, the cut-off brightness seems 
somehow to flatten over the age range 3-7 Gyr (see Fig. ll8i . 
but such plateau is incurably fainter than the observed cut-off. 
Moreover, we notice that metallicity effects are likely to affect 
the IFMR, i.e. larger final masses at lower metallicites. 

Ciardullo (2003), on the base of Ciardullo & Jacoby's 
(1999) paper, suggests that the stringet correlation between 
mass of the central star and circumstellar extinction (due to 
dusty envelopes ejected during the AGB phase) should pre- 
vent high-core mass PNe from populating brighter bins beyond 
the observed PNLF cut-off. In other words, over-luminous 
Population I PNe (evolved from young and massive progen- 
itors, hence potentially detectable as very bright [O m] /15007 
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Fig. 25. The PNLFs for 2 template galaxies: a massive elliptical 
and a spiral one. See the text for more details. 



sources in spiral and irregular galaxies) would be efficiently 
extinguished below the observed cut-off. 

To this respect we notice that there is an additional factor, 
other than extinction, that would locate high-core mass PNe at 
fainter magnitudes than the cut-off. To this respect, we should 
consider that one of the underlying assumptions in Ciardullo 
& Jacoby (1999) analysis is that a constant fraction of the stel- 
lar luminosity is reprocessed into [O m] /i5007 emission, which 
yields a positive correlation between the mass of the central 
star and the maximum [O in] /15007 luminosity emitted during 
the PN phase. Our calculations, instead, show that this correla- 
tion is not monotonic, as the maximum [O in] A5QQ1 luminosity 
presents a peak at central star masses of about 0.70 - 0.75 M , 
then declining at larger masses. This implies that, even with- 
out the effect of circumstellar extinction, PNe with high-mass 
central stars would not contribute to the PNLF cut-off because 
intrinsically fainter. 

In conclusion, we should stress that our difficulty to explain 
the age-invariance of the PNLF cut-off could not be solved by 
including the effect of circumstellar extinction, as proposed by 
Ciardullo & Jacoby (1999). In fact, we do not have any need to 
make PN models with massive progenitors fainter, as required 
by Ciardullo & Jacoby (1999), but rather to make PN mod- 
els with less massive and old progenitors several magnitudes 
brighter. How to obtain this effect with the aid of consistent PN 
models, is still to be investigated. 

5.3. The dependence on galaxy type: the age 
alternative 

The PNLF is expected to depend on the galaxy type and even 
more relevant here on the age of the stellar populations. In or- 



Fig. 26. The same as Fig.^] but now simulating SFR with just 
2 bursts of SFR: an old an intense burst in the beginning of 
the galaxy history, of age 10 Gyr, added to a more recent burst 
just 5 percent as intense as the old one. It can be seen that the 
PNLF shapes change a lot, with respect to the continuous SFR 
shown in Fig.[D3 However, the bright PNLF cutoff is found in 
the same magnitude bins as before. 

der to discuss this point, we start by presenting in Fig.l25lthe 
predicted PNLFs for two template galaxy models representing 
two extreme cases: massive ellipticals and spirals. The ellipti- 
cal galaxy corresponds to the 10 12 M model from Chiosi & 
Carraro (2002), in which the star formation is all concentrated 
in the first 2 to 4 Gyr of the galaxy life; during the same pe- 
riod its metal content rapidly increases up to near-solar values. 
The spiral galaxy model, instead, is derived from chemical evo- 
lution models of our own Galaxy: it corresponds to model B 
from Portinari et al. (1998). The star formation in this case is 
characterized by a initial rapid increase during the first 3 Gyr, 
followed by ever-decreasing, but continuous, SFR. The AMR 
is derived from models but is in excelent agreement with ob- 
servational data from disk stars. 

These template galaxy models show clearly a difference of 
about 5 mag between the cut-off position expected (and ob- 
served) in spirals at M* ^ -4.5, and the one expected in old 
ellipticals at M* ^ 0.5. This difference follows, essentially, 
from significant age dependence expected by our models (see 
Fig-CSJ, and is in clear contrast with observational claims (see 
Sect.G). 

Thus, if ellipticals were systems consisting exclusively of 
old stellar populations, we could hardly explain the invariance 
of the cut-off. A viable explanation for the observed invariance 
is that the ellipticals with measured PNLFs may have had a 
burst of star formation in recent times, with ages close to 1 Gyr. 
In order to illustrate this possibility, in Fig.[26]we present a se- 
ries of simulated PNLFs for galaxies having an intense burst 



20 



Marigo et al.: Population effects in the PNLF 



of SFR 10 Gyr ago, followed by a more recent one at the age 
fburst! which involved just 5 % of the stellar mass in the galaxy. 
It can be noticed that the shape of these PNLFs does change 
a lot, compared to the continuous-SFR PNLFs illustrated in 
Fig. ^] Actually, in the cases of bursting SFR, the PNLFs 
present a brighter component, due to the more recent burst, that 
displays a cut-off behaviour over the brightest bins, while con- 
tinuously declining towards fainter magnitudes. Then, at about 
M5007 ^ 0.3, a second, highly-populated and faint component 
builds up as a consequence of the old burst. We also notice that 
if just the brightest PNe were observed in such galaxies, the 
derived cut-offs would be exactly those caused by the youngest 
burst. Only deeper observations could then reveal the different 
shape of these brighter components, and the sharp rise - ac- 
tually a second cut-off - at M5007 ^ 0.3. Therefore, if PNLFs 
are observed deeply enough they can give very valuable results 
about the star formation history during the last several billion 
years. 

In summary, in the context of our PN models, the forma- 
tion of the PNLF cut-off, observed in some elliptical galaxies 
as bright as in spirals, could be explained only invoking a re- 
cent burst of SFR, as young as a; 1 Gyr, corresponding to stellar 
progenitors with initial mass of about 2.5 - 3.0 M , or corre- 
spondingly to CSPN mass of about 0.70 - 0.75 M . 

Let us also clarify that a change of ill-known model pa- 
rameters hardly would change this scenario. In fact, in order 
to reconcile our models with the claimed age-invariance of the 
PNLF cut-off, we would have to get much brighter PNe mod- 
els at old ages. By using He-burning post-AGB tracks instead 
of H-burning ones, our old models become even fainter than 
above mentioned. Also, older PN models are little sensitive to 
the exact value of the transition time, and choosing these times 
much larger than 500 yr - as likely appropriate for models of 
lower mass - would also cause fainter old PNe models, instead 
of brighter ones. 

The idea the early type galaxies may suffer secondary re- 
cent star formation activity has long been around in literature. 
It is hinted by the diagnostics of the line absortption indices on 
the Lick system, such as HB and/or Hy, that are commonly 
considered as good indicators of the age of the most recent 
star forming episode (Bressan et al. 1996; Jorgensen 1999; 
Kuntschner 2000; Longhetti et al. 2000; Trager et al. 2000a,b; 
Vazdekis et al. 2001; Davies et al. 2001; Poggianti et al. 2001; 
Tantalo & Chiosi 2003), Further supported by evidences from 
various broad-band colors (Zabludoff 1996; Goto et al. 2003; 
Quintero et al. 2004; Fukugita et al. 2004), and finally found 
also in NB-THSP simulations of galaxy formation and evo- 
lution (Kawata 2001a,b; Carraro & Chiosi 2003; Kawata & 
Gibson 2003a,b). However, the constraint set by the PNLF cut- 
off luminosity is very demanding because it would imply that in 
all ellipticals showing bright PNs star formation activity should 
have occurred as recently as about 1 Gyr ago in order to reach 
the same luminosity of the ongoing star formation case as in 
spirals. 

The question may be posed whether sampling ellipticals 
and spirals to derive the empirical PNLF, selection effects 
somehow favour those ellipticals with bright PNe, i.e. recent 
star formation, so that ellipticals with older stellar content and 



whence fainter cut-off luminosity for their PNLF are simply 
missed. This issue is not specifically mentioned in literature so 
that a clear answer cannot be given here. 

The most puzzling aspect with the age interpretation is 
that there are some elliptical galaxies, for instance NGC4649, 
NGC4406, NGC4486, NGC1399, and NGC1404, which have 
bright PNs but yet possess broad colors and line absorption in- 
dices - for instance mean B — V = 0.96 and mean HB - 1 .45 
- strongly indicating that only old stars are present, unless the 
last episode has involved less 0.05% of the galaxy mass as thor- 
oughly discussed by Tantalo & Chiosi (2003). 

A possible way to remove the theoretical discrepancy might 
be related to the IFMR and its dependence on metallicity. In 
practice, we could obtain much brighter PNe at old ages if their 
low-mass and metal-poor stellar progenitors had been able to 
build up relatively massive cores at the end of the AGB evo- 
lution, hence brighter CSPNe. Indeed a general trend of in- 
creasing final mass at decreasing Z is already predicted by our 
models (see Fig. I12i . but the resulting effect on the PNLF cut- 
off is not strong as the lowest metallicity here considered is 
just Z = 0.004. Additional calculations at lower metallicity are 
then required to test the above hypothesis. This analysis is in 
progress. Alternatively, we are left with the open possibilities 
that either the current models of PNs are still too simple and 
neglect some important effect (we are, however, confident that 
this is not the case), or some new type of old stars intervene 
in particular in the case of old elliptical galaxies in such a way 
they closely mimic the properties of classical PNs (investiga- 
tion of this topic is currently underway). 

6. Concluding remarks 

In Paper I we had developed a theoretical model for synthetic 
PNe. In the present paper we use this model to predict proper- 
ties of synthetic samples of PNe. In particular we investigated 
the PNLF of synthetic PN populations with regard to its shape, 
maximum [O m] A5QQ1 brightness, and how it depends on age 
and composition of the population as well as on assumptions 
necessarily entering our model. With respect to Paper I and in 
view of the importance of a reliable prediction of the emis- 
sion properties of PNe, we extended our model by including the 
photoionisation program CLOUDY, thereby eliminating one of 
the free parameters, the electron temperature T e in the nebulae. 
The detailed radiation transfer treatment confirms that the as- 
sumption of an approximately constant T e of about 10000 K in 
the dynamical part of the model is justified, at least within the 
accuracy required at this stage. Apart from free assumptions 
about age and composition of the population causing PNe, and 
about the SFH and IMF, the transition time remains the sole 
physical parameter for which we have only rough indications, 
and for which we would need additional theoretical calcula- 
tions of the superwind phase of AGB stars, and their transition 
to the post-AGB phase. We are also forced to assume spherical 
symmetry of the circumstellar matter. 

With this improved model for PNe populations we have 
simulated galactic PNe and compared the results with observa- 
tions. We concentrated on emission lines and ionisation states, 
and find excellent agreement without the necessity to vary any 
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model parameters. Based on these successful comparisons, we 
then concentrated on the luminosity function. The main results 
concerning the PNLF are: 

- Our models predict a PNLF shape towards the cut-off very 
similar to those observed in spiral and irregular galaxies. 

- In galaxies with recent star formation (either continued up 
to now, or with a burst occurred ~ 1 Gyr ago), the cut- 
off is expected at -5 < M5007 < -4, depending on model 
parameters. This is consistent with the observed values at 
about -4.5 (e.g. for the LMC). 

- In galaxies without recent star formation, the cut-off is pre- 
dicted much fainter than observed. This is due to the lack 
of stars with M > 2 M , responsible for the brightest PNe 
(M 500 7 < -3.0). 

- From the two previous points it follows that we can hardly 
explain the constancy of the cut-off observed in some el- 
lipticals (dominated by old populations) in terms of the 
claimed age-invariance. Rather, a possible way to explain 
the bright PNLF cut-off observed in some late-type galax- 
ies is just to invoke a recent burst of star formation. 

- One main parameter determining the PNLF cut-off is the 
transition time, to be constrained on the base of observed 
PN properties (e.g. 7V e , T e , ionized masses, radii, expan- 
sion velocities, etc). Our simulations indicate that, in or- 
der to recover the observed brightness of the cut-off, the 
CSPNe populating it (with masses around 0.70 - 0.75 M ) 
should have a f tr of the order of few hundred years (typi- 
cally 500 - 1000 yr). On galaxies without recent SFR, the 
cut-off position becomes little sensitive to t^. 

- The predicted dependence of the PNLF on metallicity is the 
result of a number of factors, as we need to know (i) how 
the evolutionary properties of the stellar progenitors (i.e. 
the efficiency of mass loss on the AGB, the details of the 
initial-final mass relation) vary as a function of the initial 
chemical composition, (ii) the initial main-sequence abun- 
dance of oxygen, and iii) the actual oxygen abundance in 
PNe since it may differ significantly from the initial value 
due to dredge-up episodes during the TP- AGB phase. 

Overall, our results explain in a very natural way the 
successful calibration of the PNLF using Cepheid distances 
(Feldmeier et al. 1997; Jacoby 1996). In the list given by 
Jacoby (1996), 7 out of 8 galaxies in the sample are classified 
as spiral or irregular galaxies, only NGC5253 is of E/S0 type. 
However, it is well documented that this galaxy has undergone 
a recent starbust (see, e.g. Verma et al. 2003) and thus, as for 
all the other objects star formation has taken place during the 
last billion year, and all PNLFs reach M* . In fact, the argument 
can be made even stronger and more general: in case there is a 
Cepheid distance to a galaxy, there must be a young population 
and thus the PNLF will reach the canonical cutoff brightness. 

Furthermore, our models do not support the idea of a 
population-invariant PNLF cut-off. Indeed, the suggestions for 
an age-invariant cut-off poses a serious problem to the theory. 
In order to explain it, we may need to explore the properties of 
very low-Z PNe, or to invoke a very different and difficult-to- 
imagine scenario for the evolution of PNe and their progenitor 
stars than here considered. 



Considering all these aspects, we find it worth questioning 
if the indications for the cut-off invariance are really as solid 
as claimed. Until definitive evidence is not brought, and until 
there is not a solid theoretical framework supporting such ev- 
idence, we suggest that the PNLF should still be considered 
as a potentially useful probe of stellar populations in galaxies, 
rather than a standard candle. 
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